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Abstract 

A parameter-free, mean-field, multi-orbital Hubbard model with nonspherical 
Coulomb and exchange interactions, implemented around all-electron local- 
density approximation (LDA) calculations, correctly predicts the band-gap 
energy, the absence of dimerization, and the antiferromagnetic ground state 
of halogen-bridged linear-chain Ni compounds. This approach also reproduces 
the insulating ground state and dimerization in PtX linear-chain compounds 

in agreement with experiment and previous calculations. 
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The halogen-bridged transition-metal linear-chain compounds, referred to as MX com- 
pounds because of their alternating transition-metal atoms M and halogen atoms X, form 
weakly coupled linear-chain-like structures, and are of considerable interest [|T] due to 
a rich and accessible phase diagram |2|-|T2||. More specifically, experimentalists can now 
tune the strength of the charge- density wave (CDW) and spin-density wave (SDW) J|.|7J 
and induce phase transitions from an CDW to an SDW state ||. First-principles density- 
functional calculations of neutral PtX compounds have shown clearly that (1) the mechanism 
of the dimerization and the insulating ground state in the MX chain systems is due to an 
electron-phonon coupling along the chain, and (2) the ligand structure is essential for the 
correct description of ground-state properties ||. New high-quality data on single crystals 
are challenging our understanding of the electronic structure of these low-dimensional sys- 
tems . The Ni MX compounds, which are normally in an SDW state, change to a mixed 
SDW/CDW state by a simple substitution of the counterion CI - with perchlorite CIO4 0. 
In a recent Peierls-Hubbard model study, local defect excitations in NiBr are found also to 



exhibit lattice distortion relative to the undistorted ground state [IU|. It is intersting to 
mention that in the Ni compounds the electronic structure at the vicinity of Fermi level is 
the one dimensional equivalent of the undoped cuprate superconductor materials, where the 
local spin-polarized density approximation (LSDA) calculations predicted a non-magnetic 
metallic state instead of the observed antiferromagnetic insulating state. Here again the 
LSDA failed to reproduce the ground state properties of the Ni compounds. 

In this Letter we study the electronic properties of both CDW and SDW chains by means 
of a newly developed mean-field, multi-orbital Hubbard model, called "LDA + U" , that has 
nonspherical Coulomb and exchange interactions and is implemented around all-electron 



local-density approximation (LDA) calculations with a linear muffin-tin orbital basis-set |]IS 
This model enables us to extend our previously successful LDA study of CDW MX systems 
|||9|] to also, for the first time, explain from first-principles the magnetic SDW systems. 
It predicts that the ground state for the NiX (X= CI or Br) compounds is undimerized, 
insulating, and a low-spin antiferromagnetic, in good agreement with experiment f7|. The 



large octahedral ligand-field splitting A = e g — 12 9 ~ 4 eV produces a low-spin electronic 
configuration {t2 g ) 6 {e g ) 1 of the Ni +3 . It also predicts that the ground state of the PtX 
(X=C1 or I) linear-chain compounds is insulating and dimerized in agreement with our earlier 
calculations ||. We have to stress that our calculations are based on a parameter-free model 
that uses a constrained LDA determination for the Coulomb and exchange interactions. 
All the results presented in this Letter have been obtained using the "LDA + U" method, 



which is extensively described elsewhere ||16|| . Here we stress that the main feature of the 



model is to introduce a discontinuity in the LDA one-electron potential for integer change 
in orbital occupancy. The nonsphericity of the Coulomb and exchange interactions, i.e., the 
o?-orbital occupation dependence on the orbital angular-momentum quantum number m and 
m', is important for a good description of the magnetic moment of the MX chains. 

The <i 32 2_ r 2 orbital, whose polarization dominates the magnetic moment, is also impor- 
tant in controlling the bonding, the dimerization, and the magnetism in these systems. The 
total energy and potentials include the exchange and the nonsphericity of the Coulomb d-d 
interaction: 

E = E LDA - [UN(N - l)/2 - JN(N - 2)/4] + ^ E U mm ,n ma n m ,_ a 

m,m' ,a 

~^~~<2 ^ / (JJmm 1 J mm'^)^ J mcF^"m l a (1) 
m^m' ,m' ,a 

The screened Coulomb U and exchange J parameters are calculated self-consistently 



in the supercell approximation as described in [17]. The orbital-dependent one-electron 



potential is given by the derivative of Eq. (1) with respect to the orbital occupancy n r 



Vma(r) = V LDA (r) + E^mm' ~ U eff) 



Tim/- 



+ E (Pmm' — Jmm' ~ U e ff)n ma + U e ff(- ~ (2) 

We define U e ff = U — J/2. The matrices U mm i and J mm ' are calculated from the screened 
Slater integrals F k (F°, F 2 , F 4 for d-electrons) , which are determined using the calculated 
values of U and J and the ratio of unscreened F 4 /F 2 [[y|. 



Model calculations for the MX chain compounds usually use one-dimensional one- or two- 



band systems in which only the M d 3z 2_ r 2 and the X p z orbitals are taken into account [|T2 
Experimentally, MX compounds have complicated structures of complex organic molecules 
ligands that are perpendicular to the chain axis and bonded to the M atoms. The bonding 
orbitals of the ligands are strongly hybridized with the d orbitals of the M atoms. The MX 
compounds are divided into three classes: (1) anionic, (2) cationic, and (3) neutral chains. 
The anionic chains have positive counterions so that the chains are negatively charged; in 
contrast, the cationic chains have negative counterions and are positively charged. The 
neutral chains have no counterions. We have studied a case of an anionic CDW chain 
K 4 PtI 4 PtI 6 , and three cationic CDW and SDW chains [M(chxn) 2 ][M(chxn) 2 X 2 ]X4 (chxn is 
1,2-diaminocyclohexane, C 6 Hi 4 N 2 ) with M= Ni (X= Br or CI) or Pt with X= CI. 

Cationic SDW chains [Ni(chxn) 2 ][Ni(chxn) 2 X 2 ]X4 . Depending on whether M is Pt or Ni, 
the associated MX chain shows a CDW or a SDW character, respectively [P^JTT|-|15| . For 



example, when MX is PtCl or PtBr, a Peierls instability is observed and the system is non- 
magnetic. When M is Ni and the counterion is CI or Br, these systems are antiferromagnetic 
and exhibit no Peierls instability. P,|T|,|r5|] . We approximate the real structure by a simpler 



model structure in order to make a direct electronic-structure calculation feasible. The 
simplest model structure is to replace the (chxn) 2 by four ammonia units per each M atom 
(see Figure 1). For the model crystal structure we compress the lattice along the a axis so 
that the distance between ammonia on neighboring chains are similar to those for ammonia- 
ligand compounds IJ. We emphasize that in our model all coordinates of the atoms are the 
same as in real crystal structure except for the absence of the CeHio ring. 

The local spin-density approximation (LSDA) alone can not reproduce the insulating 
magnetic ground-state of the SDW NiCl and NiBr systems. The electronic structure in 
the vicinity of Fermi level is very similar to that of the undoped cuprate superconductor 
materials, where LSDA calculations also predict a non-magnetic metallic state instead of 
the experimentally observed antiferromagnetic insulating state. This effect results from the 
LSDA magnetic transition being driven by the spin-polarization of a Stoner intraatomic 



exchange interaction I (about 1 eV), instead of the much stronger Hubbard interaction U 
(about 8 eV). 

Figure 2 shows the total density of states (DOS) and the energy dispersion along the chain 
direction TZ in the vicinity of Fermi level for the U — J — eV case. Despite the presence 
of the half-filled Ni d 3z 2_ r 2 band at the Fermi level the Peierls mechanism is not effective 
in opening the energy band gap. Instead, the small Ni-Ni distance sufficiently enhances the 
anharmonic elastic potential between Ni and CI that the electron-phonon coupling (Peierls 
mechanism) is rendered ineffective. |||. 

Figure 3 presents the Ni partial DOS at the vicinity of Fermi level decomposed into 
the various symmetries of the d channels for the calculated U and J as presented in Table 
1 together with the dimerization, the magnetization and the energy band gap for all the 
systems studied here. The LDA+U calculated ground state remains undimerized and the 
the Hubbard term in Eq. (2) opens an energy band gap of 1.8 eV, in good agreement 
with the experimental value of 1.9 eV ||15|| . Further, a low-spin antiferromagnetic state with 



a magnetic moment of 0.66 \ib is produced, which could be reduced by spin fluctuations 
neglected in the present study. The experimental value of the magnetic moment is not 
available at the present time. The low-spin antiferromagnetic state is due to the large 
octahedral ligand-field splitting A = e g — ti g » 4 eV which produces a low-spin (t2 9 ) 6 (e g ) 1 
electronic configuration of the Ni +3 ion. Figure 3 shows that, with respect to the the DOS 
projected on the Ni site, most of the weight of the {eg) 1 configuration comes from the 
d 3z 2_ r 2 orbital. Accordingly most of the contribution to the spin magnetic moment is from 
the polarized d 3z 2_ r 2 orbitals located in the vicinity of the band gap, and that only 0.05 
fiB is due to the polarization of the d x 2_ y 2 orbitals. In addition, the optical band gap of 
1.8 eV between the occupied Ni d% z 2_ r 2 and the unoccupied Pt d x 2_ y 2 is smaller than the 
antiferromagnetic gap of 2.3 eV arising from the splitting of the d 3z 2_ r 2 Pt bands. A slightly 
smaller value of U makes these two gaps comparable. 

The properties of the NiBr system are very similar to that of NiCl: an antiferromagnetic 
insulator with a smaller gap of 1.28 eV PJl5|. Our calculation produces a Coulomb parameter 



U = 6.9 eV, and an energy band gap of 1.6 eV with a spin magnetic moment value of 0.63 
1Mb- Here also the contribution to the spin magnetic moment is mainly from the Ni d 32 2_ r 2 
orbitals. 

Cationic CDW chains [Pt(chxn)2] [Pt^hxn^CyCL. For PtCl the experimental value of 
the dimerization is 4.9% (the short Pt-Cl distance 2.32 A and the long one 2.83 A) [|H}]. 

The LDA + U calculation shows the absence of an antiferromagnetic SDW solution for 
this system (see Table 1) and a dimerization of 4%. The 0.8 eV band gap in the dimerized 
state is much smaller than the experimental gap of 1.6 eV. This disagreement in the band 
gap is not surprising, since LDA + U method works best for fairly localized orbitals with 
spin-orbital ordering like NiO JD]. For delocalized orbitals, without spin polarization, like 



the Pt 5d LDA+U reduces essentially to the standard LDA. We should, in principle, use a 
model where the exchange correlation potential is non-local and dynamically screened, for 
example, the GW approximation of Hedin fl9f . 

Anionic CDW chain ^PtLjPtLj. This is a rather exotic material among the MX-chain 
family in that it does not have complicated organic ligands transversely bound to the metallic 
atom. We have chosen to study IQPt^Ptlg (KPtI for short), because of the simplicity of the 
ligands (I). The structure could be described as Ptl6 octahedra connected along the z-axis 
to form a chain, with K atoms located between the chains. While the measured structure is 



more complex |13| than the other systems computed in this Letter (including Pt-I-Pt angle 
smaller than 180° and disorder along the chain). We have modeled KPtI as a perfect linear 
chain without disorder. 

The computed Hubbard parameter U (computed as described above) is 3.6 eV, and is 
much smaller then the typical 7 to 10 eV values for the antiferromagnetic late-transition- 



metal oxides []T7|]. This small value of U allows no antiferromagnetic SDW solution to the 
LDA + U equations for KPtI. At computed 3.9% dimerization, The 0.37 eV energy-gap is 
much smaller than the experimental value of 1.0 eV. Like the PtCl case, the underestimation 
of the band gap is due to the inadequacy of LDA + U in handling delocalized electrons. 
Several differences between PtCl and KPtI arise from the fact that Pt-N distance (2.06 

6 



A) is much shorter than the Pt-I distance (2.65 A), and Pt-Pt distance in PtCl (5.16 A) is 
much shorter than the 5.91 Ain KPtl. In comparing PtCl to KPtI we found that PtCl has 
(i) a much strong 5<i-ligand hybridisation, (ii) a larger Pt d x 2_ y 2-d 3z 2_ r 2 energy separation, 
and (iii) a larger CDW reflected by the larger Pt 3 ~ 5 — > Pt 3+S charge transfer (0.29 electron 
in PtCl versus 0.14 electron in KPtI). 

We have studied examples of both anionic and cationic linear MX-chain compounds by 
means of a parameter-free, mean-field, multi-orbital Hubbard model implemented around 
the LDA. For the cationic NiX (X= Br or CI) chains, where the 3d orbitals are much more 
localized than the 5d of Pt, the predicted energy band gap is found to be very close to 
the experimental one. The model also predicts that the NiX chains have an undimerized 
antiferromagnetic ground-state in agreement with experiment. On the other hand, in PtCl 
and KPtI the band gap was found to be much smaller than the measured optical gaps. 
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percomputer time was provided by the Ohio Supercomputer and by the DOE. 
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FIGURES 

FIG. 1. Schematic structure of a single cationic M(chxn)2] [Ni(chxn)2X2]X4. The chxn 
(N2C6H14 ligand structure is represented in more detail in the insert (a), and the simplified struc- 
ture used in the model calculation is represented in the insert (b), where the cyclohexane is replaced 
by two H atoms. 

FIG. 2. LDA total DOS and energy dispersion along the chain direction TZ in the vicinity of 
Fermi level (5 k-points are used along the TZ direction). Although only one band of Ni d 3z 2_ r 2 
character crosses the Fermi level, the Peierls instability is not effective in opening an energy band 
gap. This is because the elastic anharmonic potential between Ni and CI ions along the chain is 
stronger than the electron-phonon coupling. 

FIG. 3. LDA+U DOS in the Ni site at the experimental lattice parameter decomposed into 
various d-channel symmetries: (a) for majority spin, and (b) for minority spin (The spin-up of 
a Ni atom is the down-spin of its nearest-neighbor Ni atoms and visa- versa). The dark shaded 
area represent the Pt d 3z 2_ r 2 weight, and the light shaded area the d x 2_ y 2 weight (the occupation 
numbers of the d 3 ?2_ r 2 and d v 2_.,2 orbitals are nX = 1.0, n\, = 0.39, n], = 0.65, 
and n\ = 0.60). The Hubbard term in Eq. (2) opens the energy band gap at the Fermi level 

x — y 

of 1.8 eV and produces an antiferromagnetic ground state. Almost the entire contribution to the 
spin magnetic moment of 0.66 (Ib is from d 3z 2_ r 2 orbitals, whereas only 0.05 \xb is due to the 
polarization of the d x 2_ y 2 orbitals. 
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TABLES 

TABLE I. Calculated dimerization, magnetic moment, and band gap for KPtl, PtCl, NiBr, 
and NiCl compared to the available experimental results indicated in parentheses. 

M-M (A) U (eV) J (eV) Dimerization (%) Magnetic moment (hb) Band gap (eV) 



Ptl 


(5.91 a ) 


3.6 


1.0 


3.9 (4.1 a ) 





0.4 (1.0 a ) 


PtCl 


(5.16 b ) 


3.5 


1.0 


4.0 (4.9 d ) 





0.8 (1.6°) 


NiBr 


(5.16 e ) 


6.9 


0.7 


(0 d ) 


0.63 


1.6 (1.28 e ) 


NiCl 


(4.89 e ) 


6.8 


0.7 


(0 e ) 


0.66 


1.8 (1.9 C ) 



Ref. |T| b Ref. Q c Ref. § d Ref. @ c Ref. |T 
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